function [x,y,z,E,nu]=grdecl2fem(fn,varargin)

grdecl  = scan_GRDECL(fn);


% disp(varargin)
if nargin>1
    atribfn = varargin{1};
    grdecl  = scan_GRDECL(atribfn, grdecl);

end

grdecl.COORD=round(grdecl.COORD,5);
grdecl.ZCORN=round(grdecl.ZCORN,5);
[X,Y,Z]=determCORNERPOINTcord(grdecl);

x=unique(X(:));
y=unique(Y(:));
z=unique(Z(:));

clear X Y Z

if nargin>1
    E=grdecl.YOUNGMOD;
    nu=grdecl.POISSONR;
end